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ABSTRACT 


Geometric  and  variational  techniques,  along  with  the  method  of  Lagrange 
multipliers,  optimal  control  theory,  and  elementary  calculus  are  the  tools  used  to 
derive  some  generalizations  of  Snell’s  Law.  The  focus  is  both  on  Snell’s  Law  of 
Refraction  and,  on  the  reflection  principle  under  various  conditions. 
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I.  INTRODUCTION 


Snell’s  Law  of  refraction  and  the  principle  of  reflection  have  both  been  known 
for  several  hundred  years.  The  intent  of  this  paper  is  to  develop  some 
generalizations  of  these  observed  phenomena.  We  derive  these  generalizations  by 
imposing  certain  constraints  upon  the  problem  and  show  that  in  all  cases  our  general 
results  reproduce  the  now  familiar  forms  of  the  Law  of  Refraction  and  of  the 
reflection  principle  when  suitably  restricted. 

In  this  paper  we  make  use  of  variational  methods  along  with  techniques  from 
optimal  control  theory  to  develop  a  generalization  of  Snell’s  Law  of  Refraction  in 
which  the  velocity  of  light  in  the  two  media  depends  upon  direction  alone.  We  also 
develop  several  interesting  inverse  problems  which  involve  the  determination  of 
various  indices  of  refraction  corresponding  to  certain  families  of  curves.  For 
instance,  if  we  desire  to  move  from  a  point  A  to  a  point  B  along  a  cycloid,  in  order 
to  avoid  an  obstruction  in  our  path,  we  can  determine  the  necessary  index  of 
refraction  needed  to  do  so. 

'The  Dutch  astronomer  and  mathematician  Willebrord  Snell  (1591-1626) 
discovered  the  law  of  refraction  by  considering  a  ray  of  light  that  passes  from  one 
medium,  in  which  light  travels  at  a  velocity  v„  to  a  second  medium,  in  which  light 
travels  at  a  velocity  v2."  [Ref.  1], 
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Fermat’s  Principle  states  that  the  time  elapsed  in  the  passage  of  light  between 
any  two  fixed  points  of  a  medium  is  an  extremum  with  respect  to  possible  paths 
connecting  the  points  [Ref.  2].  Since  the  distance  travelled  is  proportional  to  the 
time  within  each  medium,  the  path  of  minimum  time  in  a  medium  which  has  a 
constant  velocity  of  propagation  is  a  straight  line.  If  however,  we  are  given  two 
media  in  contact  with  propagation  velocities  vl  and  v2  (v,  *  v2),  the  path  followed 
between  a  point  (X^YJ  in  medium  1  and  a  point  (X2,Y2)  in  medium  2  consists  of  two 
line  segments  meeting  at  a  point  (Xo,Y0)  on  the  interface  of  the  two  media.  This 
"bending"  of  the  light  ray  is  called  refraction,  and  hence,  Snell’s  Law  of  Refraction, 
which  is  given  by: 

smffr,)  sin(4>2) 

V1  V2 

where  4>,  and  <P2  are  the  angles  formed  by  the  ray  and  the  normal  to  the  interface. 
The  paths  in  which  light  particles  travel  are  represented  by  rays,  which  indicate  the 
direction  of  propagation.  These  rays  are  perpendicular  to  the  wave  fronts,  regardless 
of  isotropy.  An  isotropic  medium  is  one  in  which  the  properties  of  the  medium  are 
the  same  in  all  directions  from  any  given  point.  A  medium  is  homogeneous  if  its 
properties  do  not  vary  from  point  to  point.  Wavefronts  can  have  many  shapes.  If 
the  disturbances  are  propagated  in  a  single  direction,  the  waves  are  called  plane 
waves.  "At  a  given  instant,  conditions  are  the  same  everywhere  on  any  plane 
perpendicular  to  the  direction  of  propagation."  [Ref.  l:p.  301], 
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It  has  been  observed  that  if  a  plane  light  wave  falls  on  another  plane  surface 
the  light  beam  is  reflected  from  the  surface  [Ref.  1].  The  angle  of  incidence,  and 
the  angle  of  reflection,  03,  measured  from  the  normal  to  the  surface  satisfy  the 
principle  of  reflection  which  is  given  by: 

<j>,  =  <d3 

In  the  study  of  acoustic  wave  propagation,  Snell's  Law  of  Refraction  and  the 
law  governing  the  reflection  principle  are  usually  derived  from  the  reduced 
Helmholtz  equation.  Although  this  thesis  concerns  generalizations  of  these  laws 
using  variational  methods,  a  brief  outline  of  this  process  is  included  to  illustrate 
other  available  techniques. 

Figure  1.1  depicts  plane  wave  reflection  and  refraction  from  a  planar  interface. 
The  equation  for  the  incident  plane  wave,  P,  is  given  by: 

p  =g  a.fxsmS, ->0060,; 

where  the  incident  amplitude  has  been  normalized  to  1,  the  vector  x  is  a  position 
vector,  and  3  is  the  vector  which  points  in  the  direction  of  propagation  of  the  wave. 

Using  subscripts  u  for  the  upper  medium,  1  for  the  lower  medium,  i  for  incident,  r 
for  reflected,  and  t  for  transmitted  we  have  the  following  system  of  equations: 

v:P,  +k:uP,  =  0,  y  >  0 
v2Pr  +k2uPr  =  0,  y  >  0 
v2?,  +k:;P,  =  0,  y  <  0 


Figure  1.1  Plane  Wave  Reflection  and  Refraction 

where  k  =u/c  =  wave  number.  An  e  '“'  time  dependence  is  assumed  throughout  the 
following  analysis.  Allowing  for  continuity  of  pressure  at  the  interface,  we  have  the 
condition: 

P„  =  P;  +  P,  =  P,  =  P, 
or, 
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p, +  p,  =  p, 

For  an  ideal  fluid,  we  have: 

VP=-pV4>=iG>p« 

where  u  is  the  fluid  velocity.  After  taking  the  dot  product  of  both  sides  of  this 
equation  with  n  we  have: 

n-VP=ii  -icapu 

where  n  is  the  unit  normal  to  the  interface.  This  equation  yields: 


dP  . 

—  =«*pv, 
dy 

where  v  is  the  fluid  velocity  normal  to  the  interface.  Our  second  boundary  condition 


at  the  interface  equates  normal  fluid  velocities,  yielding: 


_L_a 

Pu 


where  pu  and  p,  are  densities  of  the  two  fluid  mediums.  We  define  reflected  and 
transmitted  plane  waves  as 


and 


Pr=Re 


ikKxJr 


P^Te*1**. 

Substituting  these  into  the  two  interface  boundary  conditions  produces: 
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^  JkjjamQ,  +  ^  ikjam6r  _  Jf^rnnO. 

-cose, f^'^cose^^^'U—'cosB.Tf . 

ra  Pi 

In  order  for  these  two  equations  to  be  satisfied  for  all  x,  the  phases  must  be  equal 
along  the  interface.  From  this  we  have: 

ikusin@,  =  ik„sin0r  =  ik,sin©t. 

This  equation  yields  both  the  law  of  refraction  and  the  law  governing  reflection: 

sin0.  sin0. 


and, 

0i  =  ©, 

The  solution  to  the  problem  for  R  and  T  is  that 


r= 


2cos0. 


cosfy-aPcosO, 

cos0.+ajJcos0f 


- ,  where  a  =  —  and  P  = — -. 

cosO^  +a  Pcos0r  p,  Cj 


Special  cases  such  as  Brewster’s  angle  and  the  critical  angle  where  total  internal 
reflection  occurs  can  be  studied  in  the  context  of  how  R  and  T  vary  for  differing 
material  properties  and  indices  of  reflection. 

Reflection  and  refraction  are  both  important  topics  in  the  study  of  acoustic, 
electromagnetic,  and  elastic  wave  propagation.  The  following  list  of  books  all  have 
good  treatments  of  this  subject  matter:  Halliday  and  Resnick  [Ref.  1],  Jones  [Ref. 
3],  Lorrain  and  Carson  [Ref.  4],  and  Achenbach  [Ref.  5]. 
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II.  MATHEMATICAL  PRELIMINARIES 

A.  INTRODUCTION 

The  purpose  of  this  chapter  is  to  present  a  brief  overview  of  the  basic 
mathematical  concepts  which  are  employed  in  future  chapters.  The  important 
concepts  include  the  idea  of  a  functional,  the  simplest  variational  problem,  the 
Euler-Lagrange  necessary  condition  for  the  solution  of  the  simplest  variational 
problem,  special  cases  of  the  Euler-Lagrange  equation,  an  introduction  to  optimal 
control,  Pontryagin’s  Maximum  Principle,  and,  the  use  of  Lagrange  multipliers.  The 
following  references  were  used  in  the  development  of  this  material:  Gelfand  and 
Fomin  [Ref.  6],  Weinstock  [Ref.  2],  Leitmann  [Ref.  7],  Ewing  [Ref.  8],  Weir  [Ref.  9], 
Bliss  [Ref.  10],  Pontryagin  [Ref.  11],  Taylor  and  Mann  [Ref.  12],  and,  Fleming  and 
Rishel  [Ref.  13]. 

B.  FUNCTIONALS 

Variable  quantities  called  functionals  play  an  important  role  in  several  areas 
of  mathematics,  physics  and  mechanics.  A  functional  is  an  operator  which  maps  a 
function  onto  a  real  number. 

For  the  problem  of  wave  propagation  in  a  host  medium,  consider  all  possible 
paths  joining  two  given  points  A  and  B  in  the  plane.  "Suppose  that  a  particle  can 
move  along  any  of  these  paths,  and  let  the  particle  have  velocity  v(x,y)  at  the  point 
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(x,y).  We  obtain  a  functional  by  associating  with  each  path  the  time  required  for  the 
particle  to  traverse  the  path."  [Ref.  6:p.  1]  Theoretical  proofs  of  some  results  in 
calculus  of  variations  depend  upon  concepts  from  functional  analysis.  For 
applications  of  functional  analysis  and  variational  methods  in  engineering,  see  Reddy 
[Ref.  14].  For  a  discussion  on  applications  of  functional  analysis  to  time  optimal 
control  problems,  consult  Hermes  and  LaSalle  [Ref.  15]. 

C.  THE  SIMPLEST  VARIATIONAL  PROBLEM 

The  simplest  variational  problem  involves  the  determination  of  the  maxima  or 
minima  of  a  functional  of  the  form: 

J[y]  =  [  bF(x,y,y')dx 

with  fixed  boundary  conditions:  y(a)  =  A  and  y(b)  =  B.  An  example  of  the  simplest 
problem  is  finding  the  shortest  plane  curve  joining  two  points  A  and  B.  That  is.  the 
problem  is  to  find  the  curve,  y  =  y(x).  for  which  the  functional 

J[y]=(bs[updx 

J  a 

achieves  its  minimum.  The  curve  in  question  turns  out  to  be  a  straight  line 
connecting  the  two  points  A  and  B. 

D.  THE  EULER-LAGRANGE  NECESSARY  CONDITION 

If  there  exists  a  continuous,  twice-differentiable  function  that  minimizes  the 
simplest  problem,  what  then  is  the  differential  equation  which  satisfies  this  function? 
The  answer  to  this  question  leads  to  the  implementation  of  the  Euler-Lagrange 
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equation.  We  include  the  following  sketch  of  a  derivation  of  the  Ealer-Lagrange 
equation  for  completeness.  For  an  in-depth  development  of  this  matter,  see 
Weinstock  [Ref.  2:pp.  20-22].  We  know  that  f(x  +  h)  -  f(x)  =  f'(x)h  +  higher  order 
terms  in  h.  We  apply  this  to  our  simplest  variational  problem  by  considering  two 
functions  y  and  y  +  h  satisfying  the  same  boundary  conditions.  Let  y(a)=  A  and 
(y  +  h)(a)  =  A,  then  h(a)  =0  The  same  argument  holds  true  for  y(b)  =  B.  Now, 
define  the  following: 

aJ  =  J[y  +  h]  -  J[y] 

fyF(x,y+h,y'+h)-F(x,y,yf))dx 
and  use  a  Taylor  series  expansion  to  arrive  at: 

&J=  [  b(F  h+F Jt')dx  +  higher  order  terms.  We  also  define  <$J[h]  to  be  the  first 
J  a  *  ' 


variation: 

bJ[h\=f\Fyh+Fylh)dx. 

Setting  the  first  variation  equal  to  zero,  and  integrating  by  parts  we  obtain: 

iy ,**♦/>,*  '*-/>,**♦  v  ■  °  • 

Now,  since  h(a)  =  h(b)  =  0, 


This  leaves  us  with: 
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pfr-r'vWo- 

’  ax  _ 

Ry  the  fundamental  lemma  of  the  calculus  of  variations,  (see  Gelfand  and  Fomin 
(Ref.  6:p.  9]),  it  follows  that: 

/'  b 

a(x)h(x)dx=0,  where  h(a)  =  h(b)  =  0,  and  a  and  h  are  C[a,b]  then  a(x)  =  0 

a 

for  all  x  e  [a.b].  Then, 

F  -—IF  )  =  0. 

y  j  v y ' 

Thus,  the  Euler-Lagrange  equation  for  the  simplest  problem  is: 

o- 

Thus,  a  function  y(x)  which  satisfies  the  Euler-Lagrange  equation  and  the  associated 
boundary  conditions  is  by  definition  the  extremizing  function. 

The  use  of  the  Euler-Lagrange  equation  does  not  imply  that  the  solution  is  an 
extremum.  That  is.  a  solution  to  the  Euler  equation  presupposes  the  existence  of  a 
maxima  or  a  minima.  Just  as  in  the  ordinary  differential  calculus,  if  we  consider 
f(x)  =  x3,  setting  f ' ( x)  =  3x2  =  0  yields  neither  a  minima  or  a  maxima.  Another 
necessary  condition  is  the  Legendre  necessary  condition  which  states  that  for  the 
functional: 
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J\y]=fbF(x,y,y)dx, 

J  a 

y(a)  =  A,  y(b)  =  B  to  have  a  minimum  for  the  curve  y  =  y(x),  the  inequality 
Fyv  >  0  must  be  satisfied  at  every  point  of  the  curve.  Other  necessary  conditions 
are  further  developed  in  Leitmann  [Ref.  7],  and,  Bliss  [Ref.  10]. 

E.  SPECIAL  CASES  OF  THE  EULER-LAGRANGE  EQUATION 

The  Euler-LagTange  equation  plays  an  important  role  in  the  calculus  of 
variations  and  is  normally  a  second-order  differential  equation.  There  are  some 
special  cases  for  which  Euler’s  equation  can  be  reduced  to  a  first-order  differential 
equation.  For  instance,  if  the  functional  has  no  y  dependence  Euler’s  equation 
becomes  Fy  =  C.  If  the  functional  has  no  x  dependence  then  the  corresponding 
Euler  equation  is  F  -  y'Fy  =  C.  Finally,  if  the  functional  does  not  depend  upon  y', 
Euler’s  equation  reduces  to  Fy(x,y)  =  0  and  therefore  is  not  a  differential  equation. 
For  further  treatment  of  this  subject,  see  Gelfand  and  Fomin  [Ref.  6:pp.  18-19],  and, 
Weinstock  [Ref.  2:p.  52], 

F.  OPTIMAL  CONTROL 

The  problem  of  optimal  control  is  intimately  related  to  certain  problems  in  the 
calculus  of  variations.  The  following  is  an  outline  of  this  process  and  is  taken  from 
Gelfand  and  Fomin  [Ref.  6:pp.  218-219],  Suppose  the  state  of  a  physical  system  is 
characterized  by  n  real  numbers  x‘,...,x\  forming  a  vector  x  =  (x',...,x“)  in  the 
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n-dimensional  phase  space  X  of  the  system.  Further,  suppose  the  state  varies  with 
time  in  the  way  described  by  the  system  of  differential  equations: 


dx 1 

—  Here  the  k  real  numbers  u1 . uk  form  a  vector 

dt 

u  =  (u',...,uk)  belonging  to  some  fixed  control  region  Q,  which  we  take  to  be  a  subset 
of  k-dimensional  Euclidean  space.  The  f!(x,u)  are  n  continuous  functions  defined  for 
all  x  e  X  and  all  u  e  Q.  Now  suppose  we  specify  a  vector  function  u(t),  tc  <  t  £  t, 
called  the  control  function,  with  values  in  £2.  Then,  substituting  u  =  u(t)  in  the 
differential  equations  above  yields: 

dt 

For  every  initial  value  of  x0  =  x(t0).  this  system  has  a  definite  solution,  called  a 
trajectory.  The  aggregate,  U  =  {u(t),t0,t1,x0},  consisting  of  a  control  function  u(t).  an 
interval  [t^t,]  and  an  initial  value  x0  =  x(t0),  will  be  called  a  control  process.  Thus, 
to  every  control  process,  there  corresponds  a  trajectory.  Now  let  f°(x! . x".u'....,uk) 

be  a  function  which  is  defined,  together  with  its  partial  derivatives, 

dx‘ 


for  all  x  e  X,  and  all  u  e  Q. 

To  every  control  process  U,  we  assign  the  number  J‘l f°(x,u)dt.  Then,  the 
control  process  is  said  to  be  optima!  if  the  inequality.  J[U]  <  J[U*]  holds  for  any 


other  control  process  U*  carrying  the  given  point  x*,  into  the  point  x,.  By  the  optimal 
trajectory,  we  mean  the  one  corresponding  to  the  optimal  control  process.  The 
following  hooks  have  detailed  explanations  of  this  material  and  are  recommended 
for  further  study:  Pontryagin  [Ref.  11],  Leitmann  [Ref.  7:pp.  79-83],  and  Berkovitz 
[Ref.  16]. 


G.  PONTRYAGIN’S  MAXIMUM  PRINCIPLE 

We  now'  state  the  main  theorem  of  optimal  control  theory,  the  maximum 
principle.  For  further  reading  on  this  theorem  and  its  proof,  see  Pontryagin  [Ref. 
ll:p.  11],  Gelfand  and  Fomin  [Ref.  6:p.  222],  and  Leitmann  [Ref.  7:p.  118]. 

Let  U  =  {u(t),tu,t„x0}  be  an  admissible  control  process,  and  let  x(t)  be  the 
corresponding  integral  curve  of  the  system, 

dx 1 

- =f‘(x,u)( i  =  0,1 . n),  passing  through  the  point  (O.xV . x^)  for  t  =  0,  and 


satisfying  the  conditions:  x'(t,)  =  x;' . xn( tt)  =  x,n  for  t  =  t,.  Then  if  the  control 

process  U  is  optimal  there  exists  a  continuous  vector  function  q<(t)  =  (t0(t), 
*,(t) . 'Mt))  such  that: 


1.  The  function  >t(t)  satisfies  the  system, 
(i  =  0.1 . n).  for  x  =  x(t).  u  =  u(t). 


dV. 

dt 


n 

-E 

a=o 


dT(x,u )  y 

ax' 
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n 


2.  For  all  t  in  [t^t,],  the  function  Y^Oc,!*) achieves  its  maximum 

a=0 

for  u  =  u(t). 


3.  The  relations,  #0(ti)  s  0,  G[4'(t,),u(t!)]  =  0  hold  at  the  time  t„  where 

g(*  ,!)=*“£  Her  m). 


H.  LAGRANGE  MULTIPLIERS 

"A  necessary  condition  for  a  minimum  (or  maximum)  of  F(x,y . z)  with  respect 

to  variables  x,y,...,z  that  satisfy  G,(x,y,...,z)  =  C„  (i  =  1,2 . N),  where  C  are  given 

constants  is: 

dF*  _  dF •  _  _  dF' 
dx  dy  dz 

N 

where  F*=F+£  The  constants  a„  a2 . as,  introduced  as  undetermined 

C-l 


Lagrange  multipliers,  are  evaluated  together  with  the  minimizing  (or  maximizing) 
values  of  x,v,...,z  by  means  of  the  set  of  equations  consisting  of  those  listed  above.” 
[Ref.  2:p.  6], 

Taylor  and  Mann  [Ref.  12],  and  Weir  [Ref.  9]  give  good  treatments  of  the 
Lagrange  multiplier  method  with  several  examples. 
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III.  REFLECTION 


A.  INTRODUCTION 

Mathematical  problems  are  frequently  inspired  through  observations  of  various 
physical  phenomena.  In  fact  several  solutions  to  mathematical  problems  are  inspired 
by  nature.  The  physics  of  the  problem  provide  us  the  clues,  without  which  we  would 
have  very'  little  chance  of  succeeding  in  finding  solutions.  One  of  those  problems 
suggested  by  nature  involves  the  reflection  of  light  particles  on  a  plane  mirror. 

We  know  from  Fermat’s  Principle  that  a  light  particle  will  choose  the  path 
which  will  get  it  to  its  destination  in  the  shortest  amount  of  time.  We  also  know  that 
the  shortest  distance  between  two  points  is  a  straight  line.  But  what  happens  to 
these  light  particles  when  they  undergo  a  reflection?  Do  they  still  choose  the 
shortest  path  and  if  so,  what  is  the  shortest  path?  This  chapter  answers  these 
questions  and  others.  The  material  in  this  chapter,  although  well  known,  is 
presented  for  completeness.  We  review  the  classical  results  of  reflection  and 
complete  our  discussion  with  a  derivation  of  the  reflection  principle  using  the 
method  of  Lagrange  multipliers.  The  following  references  were  used  extensively  in 
the  preparation  of  this  chapter:  Courant  and  Robbins  [Ref.  17],  Polya  [Ref.  18], 
Kazarinoff  [Ref.  19],  Williams  [Ref.  20],  and  Chakerian  and  Ghandehari  [Ref.  21]. 
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B.  HERON’S  PRINCIPLE 


Our  problem,  as  shown  in  Figure  3.1  is  as  follows:  given  two  points  A  and  B, 
both  on  the  same  side  of  a  straight  line  L,  all  in  the  same  plane,  what  is  the  point 
C  on  L  which  minimizes  the  sum  of  the  distances  from  A  and  B? 


Figure  3.1  Heron’s  Problem 

That  is,  find  the  point  C  such  that:  AC  +  BC  is  a  minimum.  If  we  can  show  that: 

AC'  +  BC'  >  AC  +  BC 
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then  we  prove  that  the  point  C  does  in  fact  minimize  the  sum  of  the  distances  from 
A  and  B  to  C.  We  know  that: 


AC'  +  BC'  =  AC'  +  B'C' 

for  B'  the  point  reflected  in  the  line  L  through  the  perpendicular  to  L  at  E.  This 
follows  since  the  triangles  BEC'  and  B'EC'  are  congruent  triangles  and  the  line  L 
is  the  perpendicular  bisector  of  the  line  segment  BB'.  Now  by  the  triangle 
inequality,  we  have: 

AC'  +  B'C'  >  AB' 

Looking  at  Figure  3.1,  we  can  see  that: 

AB'  =  AC  +  CB' 

Now  by  the  same  argument  used  above  (congruent  triangles),  we  know  that: 

CB'  =  CB. 


Then  it  follows  that  : 


AB'  =  AC  +  CB. 


Therefore, 

AC'  +  BC'  >  AC  +  BC 

as  we  desired.  So  we  have  in  fact  established  that  the  point  C  does  minimize  the 
sum  of  the  distances  from  A  and  B.  It  follows  from  geometry  that  the  angles  <t>,  and 
are  equal  and  thus  w-e  have  the  law  which  governs  the  principle  of  reflection: 

cb  =  <b, 

Therefore,  the  reflected  ray  of  light  takes  the  shortest  possible  course  between 
the  two  points  A  and  B.  This  reflected  ray  stays  within  the  plane  of  incidence  since 
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any  other  possible  direction  of  reflection  would  not  produce  the  shortest  possible 
path  from  the  point  A  to  the  point  B.  This  discovery  is  due  to  Heron,  the 
Alexandrian  scientist  of  the  first  century  A.D.  [Ref.  17].  "Simple  as  it  is.  Heron's 
discovery  deserves  a  place  in  the  history  of  science.  It  is  the  first  example  of  the  use 
of  a  minimum  principle  describing  a  physical  phenomenon."  [Ref.  18] 

Heron’s  reflection  principle  can  be  generalized  in  Minkowski  spaces.  A 
Minkowski  space  is  simply  a  finite  dimensional  normed  linear  space.  A 
generalization  of  Heron’s  reflection  principle  using  Lagrange  multipliers  in  a 
Minkowski  plane  and  related  applications  to  conics  is  given  in  Ghandehari  [Ref.  22). 
For  a  generalization  of  Steiner’s  problem  (Fermat’s  problem)  see  Chakerian  and 
Ghandehari  [Ref.  21]. 

Consider  the  following  problem:  given  an  area  A  and  one  side  of  a  triangle,  c, 
which  triangle  is  the  one  for  which  the  sum  of  the  other  sides  a  and  b  is  the 
smallest?  See  Figure  3.2.  [Ref.  17:p.  332] 

Prescribing  the  side  c  and  the  area  A  of  a  triangle  is  equivalent  to  prescribing 
the  side  c  and  the  height  h  on  c,  since  A  =  y2(ch).  Therefore,  our  problem  is 
reduced  to  finding  the  point  R  such  that  the  distance  from  R  to  the  line  PQ  is  equal 
to  the  height,  h  and  such  that  the  sum  a  +  b  is  a  minimum.  From  the  first  condition, 
it  follows  that  R  must  lie  on  a  line  parallel  to  the  line  PQ  at  a  distance  h.  This  R 
is  given  by  Heron's  Theorem  for  the  special  case  where  P  and  Q  are  equally  distant 
from  L. 
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R'  R 


Figure  3.2  Application  of  Heron’s  Principle 


C.  STEINER’S  PROBLEM 

The  following  example  comes  from  Polya  [Ref.  18:p.  145],  and  illustrates  the 
use  of  the  reflection  principle.  It  is  commonly  known  as  Steiner’s  problem  named 
after  Jacob  Steiner,  the  famous  representative  of  geometry  at  the  University  of 
Berlin  in  the  early  nineteenth  century.  Three  towns  intend  to  construct  three  roads 
to  a  common  traffic  center  which  should  be  chosen  so  that  the  total  cost  of  the  road 
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construction  is  minimized,  assuming  cost  is  proportional  to  length  alone.  Simply  put, 
given  three  points,  find  a  fourth  such  that  the  sum  of  the  distances  from  the  three 
points  to  the  fourth  is  minimized.  See  Figure  3.3.  Let  A,  B,  and  C  denote  the  three 
towns  and  X  the  fourth  point  in  question.  Also,  assume  that  all  interior  angles  a,  b, 
and  c  are  all  less  than  120°. 


Figure  3.3  Steiner’s  Problem 
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We  want  to  minimize  the  following  expression:  C  =  AX  +  BX  +  CX.  If  we  fix  one 
of  these  distances,  say  CX,  then  we  have  a  problem  which  resembles  Heron’s 
problem.  Therefore,  we  let  CX  =  r  and  minimize  AX  +  BX.  In  this  case,  we  are 
dealing  with  reflection  in  a  circular  mirror.  We  know  that  light  particles  move  in 
such  a  way  that  the  angle  of  incidence  is  equal  to  the  angle  of  reflection.  So  in  our 
problem  the  angle  <  AXB  must  be  bisected  by  a  straight  line  passing  through  C  and 
X.  Similarly  the  angles  <AXC  and  <BCX  must  also  be  bisected  by  lines  through 
B  and  A  respectively.  As  such,  the  three  straight  lines  joining  X  to  A  ,  B,  and  C 
dissect  the  plane  into  six  angles,  the  common  vertex  being  X.  It  then  follows  that 
all  six  angles  are  equal  and  equal  to  60°.  Therefore  the  three  roads  diverging  from 
the  traffic  center  do  so  at  120°  angles.  If  one  of  the  interior  angles  is  at  least  120°, 
then  the  common  vertex  which  minimizes  the  sum  of  the  distances  from  the  three 
points  A,  B,  C  will  be  that  vertex  which  has  an  interior  angle  of  at  least  120°. 

D.  REFLECTION  PROPERTIES  OF  CONICS 

1.  Ellipse 

The  reflection  principle  applied  to  an  ellipse  indicates  that  the  focal  radii 
make  equal  angles  with  the  tangent  line.  The  following  derivation  is  taken  from 
Kazarinoff  [Ref.  19].  Let  L  be  a  tangent  to  an  ellipse  with  foci  F,  and  F2,  see  Figure 
3.4.  Let  P  be  the  point  of  tangency,  and  let  Q  be  any  other  point  on  L.  Since  Q  is 
outside  of  the  ellipse: 
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QF,  +  QF2  >  PF,  +  PF2. 

Therefore,  F,PF2  is  the  shortest  path  from  F,  to  L  to  F2  and  the  angles  a  and  B  are 
the  same. 


Q 


Figure  3.4  Reflection  Properties  of  an  Ellipse 
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2.  Parabola 


An  important  property  of  parabolas  is  that  any  ray  of  light  emanating  from 
the  focus  of  a  parabolic  mirror  is  always  reflected  along  a  line  parallel  to  the  axis 
of  the  parabola.  This  property  is  equivalent  to  saying  that  the  angles  a  and  fi,  see 
Figure  3.5,  formed  by  lines  through  the  focus  F  and  an  arbitrary  point  P  on  the 
parabola  and  through  P  parallel  to  the  axis  of  the  parabola  are  equal.  For  a  proof 
of  this  property,  refer  to  Williams  [Ref.  20]. 


Figure  3.5  Reflection  Properties  of  a  Parabola 
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3.  Hyperbola 

The  tangent  properties  of  the  ellipse  are  related  to  those  of  the  hyperbola. 
Figure  3.6  depicts  a  hyperbola  with  P  and  Q  as  its  foci.  The  tangent  line  L  passing 
through  an  arbitrary  point  R  on  the  hyperbola  bisects  the  angle  subtended  at  that 
point  by  the  foci  of  the  hyperbola.  [Ref.  17:p.  335]  Therefore,  the  angles  a  and  B 
are  the  same. 


Figure  3.6  Reflection  Properties  of  a  Hyperbola 
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E.  A  GENERALIZATION  OF  REFLECTION 


We  now  turn  our  attention  to  a  derivation  of  the  law  governing  the  principle 
of  reflection  involving  a  ray  of  light  passing  through  a  homogeneous  medium,  see 
Figure  3.7.  We  make  use  of  the  method  of  Lagrange  multipliers  in  developing  our 
solution. 


Figure  3.7  Reflection  in  a  Homogeneous  Isotropic  Medium 


We  know  that: 


cos4>j 


OA  = — - —  =  asectfr. 

COs4>j 

Similarly,  for  <!>,: 

OB  =  bsec(<J>3) 

Applying  Fermat's  Principle,  we  wish  to  minimize  the  following  expression: 
asec(d\)  +  bsec(<J>3),  subject  to  the  constraint  that: 
atan(O0  +  btan(<b3)  =  c. 

Now  construct  a  new'  function,  call  it  h  such  that: 


A(< j>,,<|>3,X)  =asec(<^1)  +ftsec(4>3)  -  A.[fltan(<j>,)  +ftan(<t>3)  -c] 
Applying  the  method  of  Lagrange  multipliers  [Ref.  9],  we  set: 

dh  dh  _  dh  _q 
3<)>a  84>3  dk 

which  yields  the  following  system  of  equations: 

dh 

- -  asec^ijtanfJ),  -  Xasec^j  =  0 

34>j 

=  teec<t)3tan4>3  -  =0 

dh 

—  =  - -  irt an<t»3  -  c  =  0 
3A. 

which  yields: 
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We  then  have: 


asin<t>,  _  Xa 

COS2<j)1  COS^J 

tein4>3  \b 

COS24>3  cos2^ 

-otan<j>,  -6tan4>3  -c  =0 

sin4>j  =X  and  sin<t>3=X 


which  yields  the  law  governing  the  principle  of  reflection: 


4>!  =  4>j- 
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IV.  REFRACTION 


A.  INTRODUCTION 

"A  stick  bobbing  in  the  water  looks  sharply  bent.  We  reason  that  the  light  that 
follows  a  straight  course  in  the  water  as  in  the  air  undergoes  an  abrupt  change  of 
direction  in  emerging  from  the  water  into  the  air."  [Ref.  18]  This  is  the 
phenomenon  of  refraction.  The  light  that  travels  with  a  known  velocity  through  air, 
travels  with  a  different  velocity  through  water.  Such  a  difference  in  velocity  explains 
in  part  the  phenomenon  of  refraction.  This  chapter  will  deal  explicitly  with  Snell’s 
Law  of  Refraction  and  is  the  focus  of  this  paper.  We  develop  several  derivations, 
the  most  important  of  which  is  one  in  which  the  velocity  of  light  depends  upon  the 
direction  the  light  takes  in  each  medium.  We  use  the  method  of  Lagrange 
multipliers  and  optimal  control  theory  to  derive  this  generalization.  The  following 
references  were  useful  in  the  preparation  of  this  chapter,  Apostol  [Ref.  23],  Polya 
[Ref.  18],  Weinstock  [Ref.  2],  Taylor  and  Mann  [Ref.  12],  Pontryagin  [Ref.  11], 
Gelfand  and  Fomin  [Ref.  6],  and,  Halliday  and  Resnick  [Ref.  1]. 
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B.  A  SIMPLE  EXAMPLE 


Our  first  development,  as  an  introduction,  comes  from  A.  V.  Baez,  which  was 
published  by  Apostol  [Ref.  23]. 

It  is  known  that  of  all  possible  paths  connecting  two  points  A  and  B  in  space, 
a  light  ray  leaving  A  chooses  the  path  which  enables  it  to  reach  B  in  the  least  time. 
This  is  Fermat's  Principle  from  which  the  law  of  refraction  and  the  reflection 
principle  follow. 

Consider  the  following  problem:  A  man  is  in  a  boat  at  a  point  P,  one  mile  from 
the  nearest  point  A  on  the  shore.  He  wishes  to  go  to  a  point  B  which  is  farther 
down  the  shore,  M  miles  from  A.  He  also  knows  that  he  can  walk  faster  than  he 
can  row  the  boat.  If  he  can  row  r  miles  an  hour  and  walk  w  miles  an  hour,  toward 
what  point  C  should  he  row'  in  order  to  reach  B  in  the  least  amount  of  time?  See 
Figure  4.1.  Looking  at  Figure  4.1  if  we  let  AC  =  x,  Snell's  Law  of  Refraction 
indicates  that: 

sin<t>,  r 
sin4>2  w 

Let  <b:  =  90°  and  use  geometry  to  obtain: 

x  r 

s[U?  w 

Solving  for  x  we  have  our  solution: 

^  r 

\Jw2-r2 
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Figure  4.1  Application  of  Law  of  Refraction 

This  example  serves  to  illustrate  the  fact  that  we  can  analyze  and  solve  certain 
problems  avoiding  the  use  of  calculus  by  simply  applying  Snell’s  Law  of  Refraction. 
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C.  GENERALIZATIONS  OF  SNELL’S  LAW  OF  REFRACTION 


1.  Optically  Isotropic  Media 

Our  first  derivation  of  Snell’s  Law  of  Refraction  involves  the  passage  of 
light  between  two  points  in  two  optically  isotropic  media.  See  Figure  4.2.  That  is, 
the  velocities  of  the  particles  in  both  media  are  constant  but  different  [Ref.  2]. 
Consider  for  example  the  path  of  a  light  particle  emanating  from  a  source  at  (X„Y,) 
in  water  through  a  thin  plate  of  glass,  to  a  point  (X2,Y2)  also  in  wat°r. 


Figure  4.2  Refracted,  Reflected,  and  Incident  Rays 
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Applying  Fermat’s  Principle  from  the  point  (X„Y,)  in  medium  1  to  a  point 


(X2,Y2)  in  medium  2,  with  respective  light  velocities  of  \\  and  v2  we  arrive  at: 

T_v'(V*o>Mr,->,o)2  , 

V.  V2 

where  T  =  time  =  arc  length/velocity  =  ^ 

~dt 


Again  Fermat’s  Principle  claims  that  the  light  particle  will  take  the  path  of  shortest 
time  between  two  points  or. 


dT 


X0~X2 


V* 0 


dX, 


o  v2)jix0-X2)2  +(Y0-Y2)1  v^-Xf+^-Yj1 


X,  -xa  xn-x. 

Since  sin<t>,  =  —  -  jind  simfr.,  = 


Then  the  expression  becomes: 


sin4>j  sin4>2 


V1  V2 


2,  Optically  Isotropic  Media  Using  Lagrange  Multipliers 

We  now  attempt  to  solve  the  same  problem,  that  is  the  propagation  of 
light  in  two  isotropic  media  using  the  method  of  Lagrange  multipliers  [Ref.  12].  See 
Figure  4.3. 


Looking  at  Figure  4.3  it  is  easy  to  see  that  cos(4>,)  =  -^-  or, 

OA 


OA  =  asec  (<£,)•  Now  applying  Fermat’s  Principle  we  want  to  sum  the  arc  lengths, 

divided  by  their  respective  velocities,  to  arrive  at  an  expression  for  time: 

asec(J>j  &sec4>2 
V,  V2 

subject  to  the  constraint  that: 

atan(d>,)  +  btan(«J>2)  -  C  =  0 


Figure  4.3  Propagation  of  Light  in  Two  Isotropic  Homogeneous  Media 
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Call  the  constraining  function  giO,.^),  that  is  let: 

g(<J>!,<l>2)  =  atan(4>;)  +  btan(^>2)  -  C  =  0 
By  the  method  of  Lagrange  multipliers  [Ref.  9],  we  form  a  new  function  h  where 
h  =  T  -  AgfO),,^) 

Setting  =-^-  =  —  =0  yields  the  following  three  equations: 

3(1),  34>2  3A. 
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3.  Optically  Anisotropic  Media  With  Velocity  a  Function  of  Depth 

Now,  with  the  help  of  the  calculus  of  variations  we  consider  the  case 
where  the  velocity  is  not  constant  but  rather  a  function  of  depth.  In  other  words, 
consider  the  propagation  of  light  in  two  optically  anisotropic  media  where  the 
velocity  of  the  light  particles  is  some  function  of  y  [Ref.  24]  and  [Ref.  26],  Again, 
using  Fermat’s  Principle,  we  arrive  at  the  following  equation  which  we  want  to 
minimize: 


Jxt  ds  JX[  v(y) 

dt 

This  functional  is  of  the  form  F(y,y'),  therefore,  it  represents  a  special  case  of  the 
Euler-Lagrange  equation  when  the  integrand  has  no  x  dependence,  [Ref.  6]  and  [Ref. 
11].  Starting  with  the  basic  Euler-Lagrange  equation: 


We  derive  an  expression  which  will  simplify  this  second  order  differential  equation 


into  a  first  order  differential  equation: 


(F'}.Fy-Ft,>y'-F!Vy 


II 


Then  multiplying  by  y': 


F  v' -F  ,  \a  -F  ,  v/v//  =  -^-TF-v/F  ,1 

>y  >yy  >>'y  dx  y 

Recalling  that  a  first  integral  of  a  system  of  differential  equations  is  a  function  which 
has  constant  value  along  the  solution  curves  of  the  system,  we  can  then  say  that  the 
first  integral  of  Euler’s  equation  is: 

F  -  y'Fy'  =  C 


Using  this  equation,  we  now  solve  the  differential  equation: 

y'y'  _r 


v(y) 

From  Figure  4.2,  y^tana,  where  <j>. 


(v(y)) 

=  —  -a.  Hence, 
2 


or. 


therefore. 


sec2  a _ tan2  a 

seca(v(y))  seca(v(y)) 


COStt  £ 
v 


cos(^-  -4>j) 

- - - =C 

v 


from  which  we  arrive  at  our  familiar  expression  of  Snell’s  Law.  mainly  that: 
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sm(<j>,)  sin(4>2) 


v 


1 


4.  Optically  Anisotropic  Media  With  Velocity  a  Function  of  Direction 

Consider  the  propagation  of  light  in  two  anisotropic  media  where  the 
velocity  of  the  light  particles  is  some  function  of  direction.  That  is  v,  =  v;(a)  and 
v,  =  v2(8),  where  a  and  8  are  the  angles  formed  by  the  incident  ray  and  the 
horizontal,  and  the  refracted  ray  and  the  horizontal,  respectively,  see  Figure  4.4.  We 
solve  this  problem  using  variational  methods  and  show  that  the  result  is: 


1 


vi(a) 


-3— Ivj(a)sin(a)]  =  -^[v2(p)sin(p)]  =  C 

da  v22(p)  dP 


which,  for  the  simple  case  of  constant  velocities  (homogeneous  media)  reduces  easily 
to  our  standard  form  of  Snell’s  Law: 

sin<j>,  sin<t>2 


V1  v2 

Let  \"  =  tan(a),  then  a  -  tan  ’(y').  Applying  Fermat's  Principle,  we  wish  to 


minimize: 
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Figure  4.4  Propagation  of  Light  in  Two  Anisotropic  Media 
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7=  f iHZ-dx. 

Jxi  v,tan',(y/) 

Let - - - =/(y/).  This  is  done  co  simplify  the  arithmetic  involved.  Later,  we 

v,(tan'1Cy/)) 

replace  f(\")  by  - - - . 

v1(tan',(y/)) 

Our  problem  is  now  to  minimize: 

T=  NiT7/(yV*. 

J*1 

This  functional  is  of  the  form  F(y.y'),  therefore  we  again  make  use  of  the  special 
case  of  the  Euler  -  Lagrange  equation,  [Ref.  6]  and  [Ref.  11]. 

Euler’s  equation  then  is: 

F-y'Fy  =  Cor, 


fty,*l\+ya-y'\f(y,)>i\+y'l  + 


—^—Wy))]=C 


which  further  reduces  to: 


V1  *ya 

This  by  itself  is  a  more  general  expression  than  that  obtained  by  the  method  of 
Lagrange  multipliers  [Ref.  26].  Replacing  our  earlier  definition  of  f(y')  and  y'  we 


now  have: 


- *  —  -  tan  a  /Tuan2  a  (-—( — - — ))  =  C 

Vj(a)yT+tan2a  ^ ^ '  vi^a^ 

— VW(«)i 

1  I 


cos  a 

Vj(a) 


tana  sec  a 


v?(a) 


=  C 


or. 


cos  a 

^00 


sing 
cos2  a 


-cos2a(v,/(a)) 

Vj(a) 


=  C. 


Now  establishing  a  common  denominator  of  v;J(ar): 

Vj(a)cosa  +sina  Vj(a) 

vi(a) 

from  which  it  can  be  seen  that'. 


-iv. 

.2—  da 1  1 


Vj(a) 


(a)sin«1j 


or. 


— ! — — 
v,2(a)  *■ 


v,(a)sina]  =  — — — [v2(p)sinp]  = 
J  v22(p)rfPl 


=  C. 


We  now  desire  to  reach  the  same  generalization  of  Snell's  Law  of 


Refraction  using  optimal  control  theory. 

We  first  include  the  derivation  of  the  general  problem,  using  optimal 
control.  For  a  further  explanation  of  this  process,  see  Pontryagin  [Ref.  11].  It  is 
required  to  join  two  points  (x;,y.)  and  (x:,y2).  see  Figure  4.4,  with  a  continuous, 
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H  =  #„F(x,y,y')  +  +  *2y ' 

To  maximize  H.  we  set  —  =0: 

du 


dH  .  dF  ,  n  .  / 

— =  Ro^r +R2  =  °.  « =  y 

dii  du 


♦2  =  -Re¬ 


setting  H  =  0  we  have: 


0  =  ^0F(x,>,«)  +  tjr  ,  -(♦,—-)« 

du 

Ri  “  ~^oF(x,y,u)  +■  i|i0( — -)y' 

dy 

From  the  jump  conditions,  **  =  #'  +  ,«  grad  g  (x,y): 


♦„  -fs—fyi  - *0  -h*—fy-y 

dy  J  [  dy 


df7  df. 

-f2  *  ~(y  V  =  -/1 + -—(y  ')' + iiJV7 
dy'  dy' 


d/1  df, 


where  (N!,N:)  is  a  normal  vector  to  the  curve  g(x,y)  =  0  at  the  break  point  of  the 
trajectory.  Let  Y'  denote  the  slope  of  the  tangent  to  the  curve  of  g(x,y)  at  that 
point: 
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from  which: 


j*  _ 

_ dy'  dy' _ 

s>-A'%<yy-%<yy 


y'=-/2v1^(y)'^(ry 

dy  dy  dy  dy 


dy'  cy' 


Since  the  velocity  of  the  light  particles  is  a  function  of  direction,  the  corresponding 

integral  is:  fXif(y/)  \Jl+y'2  dx,  where /(y7)  = - — — . 

Jx >  Vjtan’fy') 


Applying  the  results  derived  above,  we  have  the  following  conditions: 

dF  n  , 

- =  —  =0  and 

dt  dx 

atj,2  _  dF  =  Q 
dt  dy 

Therefore,  4',  =  C„  and  #2  =  C2.  Using  this  information,  we  now  have: 


*  .».c 

*2  dy'  2 
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_  dF  ,  r 

'th-F'—y  =<V 

dy' 


By  the  continuity  of  we  have: 


/2~ 0'*)'=/.~0’‘)'- 
dy  dy 


We  rewrite  the  right  hand  side  of  this  equation  as: 


or. 


4^-y^yrv. 

/I +>,/2 

Recalling  that  y'  =  tan  a  and  that: 


av, 

^l  =  _5L 

3y'  v,(a)2 

we  then  arrive  at  the  expression: 


-cos2  a— 

da 

Vj(a)2 


cos  a 
v,(a) 


cos^— 

- — tanaseca 

v/a)2 


which  then  becomes: 
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v.(o)cosa  +  —  sina 

1 _ da _ 

‘.(a)2 

This  is  our  familiar  expression  for  this  particular  generalization  of  Snell's  Law  of 
Refraction: 


1  d  [v,(«)sin«]  .  — L-^jv^sinpi. 


Vj(a)2  da 


v2(P)2  dp 


5.  Optically  Anisotropic  Media  Without  the  use  of  Lagrange  Multipliers 

Consider  again  the  same  case  where  the  velocity  of  the  light  particles 
depends  upon  the  direction  at  each  point.  We  wish  to  reach  a  generalization  of 
Snell's  Law  without  the  use  of  Lagrange  multipliers.  A  look  at  Figure  4.3  shows  that 
we  need  to  minimize  the  following  expression: 

osec<t>j  bsec4>2 

where  r  =  f(<L)  is  the  polar  equation  of  a  smooth  starshaped  curve  and  r  =  g(d>;) 
is  similarly  defined.  A  smooth  starshaped  curve  has  the  property  that  every  line 
segment  connecting  any  point  within  the  curve  to  the  center  of  the  curve  remains 
completely  within  the  curve.  Another  look  at  Figure  4.3  reveals  that 

C  =  atan(<$.)  +  btan(d>2). 

Let  h(<t>..4>2)  =  atan(d>.)  +  btan(0;)  -  C.  Recalling  the  definition  of  the  total 
differential  for  a  function,  h.  where  h  =  h(d>  we  have: 
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<f4>, 

a*,  '  a»2  2 

Then  the  following  relationship  holds: 

-dk 

<*4>2  _  _  -asec2^) 

^4>!  —  isec2(4>2) 

84ij 

Now  since  we  wish  to  determine  the  minimum  of  a  function  we  simply  set 
the  derivative  of  the  function  equal  to  zero,  that  is: 


This,  and  the  chain  rule  give  us: 


dT 

dfr, 


=  0. 


asec^jtan^/^i)  -  ■taectfr/'Qfe,)] 
/(<>>!> 


feec<J)2tan<|)2g(<()2)^  ~bst ctyg'ftj 


<*2 


£2(4>2) 


=  0 


which  further  reduces  to: 


which  yields: 


asec^jtaiuljj/fvj)  -flsec<|)/'(4>,) 


asec  4), 


A*,) 

sin^zgC^)  -  cos4>2g/(«|>2)| 
8\4>2> 


/(<{>)  )sin<|>1  -/y(4>1)cos4)1]  [g(tt>2)sin4>2  -  g  /(<j>2)cos<t>2 

/($,)  £2(<t>2) 

This  expression  is  yet  another  representation  of  Snell's  Law  of  Refraction  for  our 


given  problem. 
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6.  Optically  Anisotropic  Media  in  Two  Dimensions 


Let  us  now  consider  the  problem  of  the  propagation  of  light  in  a  two 
dimensional  anisotropic  media.  That  is,  consider  a  light  particle  traveling  from  a 
point  A  to  a  point  B.  Let  the  velocities  of  the  light  particles  be  some  function  of 
two  spatial  variables,  say  y  and  z  :  v  =  v(y,z).  For  instance,  suppose  we  arrange  four 
glass  containers,  two  on  top  of  each  other,  and  place  different  liquids  in  each 
container,  resulting  in  varying  particle  velocities  in  each  medium.  See  Figure  4.5. 


Figure  4.5  Propagation  of  Light  in  Two  Dimensions 
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Applying  Fermat's  Principle,  we  arrive  at  our  familiar  expression  using  variational 
methods: 

r=  =  (Bfi*ya+zadx 
Ja  Ja  v(yj.) 

dt 

Our  functional 

J\y ]  -  [BP(yy*x/,ydx 

J  A 

has  no  x  dependence.  Therefore,  starting  with  the  two  Euler  -  Lagrange  equations 
for  our  situation: 

F,~(F)=0 
y  dx  7 

and 

F,--f(Fz,)-0. 

dx 

We  again  reduce  these  two  second  order  differential  equations  as  before 
and  are  left  with: 

F/-Fy,yya-Fyyy'y"  -  £(F-y'Fy) . 

Thus,  in  this  case.  Euler’s  equation  has  the  first  integral: 

F  -  y'Fy.  =  C 

where  C  is  a  constant.  Similarly  for  Z  we  arrive  at: 

F  -  z'F,  =  C. 

Applying  these  two  Euler-Lagrange  equations  to  our  functional  we  see  that: 
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and, 


y/l+y^+z'2  yt  y' 

v0'»z)  \/l  +ya  +Z72  v(y,z) 


y/l  +ya  +  za  _./ _ ^ _ 

v(^)  *  1  +y/2  +z^  v^) 


»C, 


which  yield: 


and. 


1  +z 


a 


V^l  +y/2+z/2  v(y,z) 


=  C, 


i+y 


\/T+/*+z'*  v(y,z) 


=  C2. 


Dividing  the  above  two  expressions  and  letting  —  =  C3w 

C2 


or. 


1  +/* 


1  +  z'2  =  C3  (  1  +  y'2). 


If  we  let  z'  =  tan(a)  and  let  y'  =  tan(B)  where  4>j  =  -<xand^2 


1  +  tan2(a)  =  C,  (1  +  tan2  (B)) 


or. 


sec2! a)  =  Cjsec2  (B) 


arrive  at 
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or  letting  C4  =  JC, 


which  in  turn  equates  to: 


sec  (a)  =  C,  sec  (8) 


cos  a  =  C.  cosp  where  C.  -  —  . 


Remembering  that  <j>j  =  —  -  a  and  <j>2  =  —  -  Pwe  now  have  the  familiar  form  of 

2  2 


Snell's  Law: 


sin(<t>,)  =  C5  sin(02). 


where  C5  represents  the  ratio  of  the  velocities  in  the  different  media. 
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V.  APPLICATIONS 


A.  INTRODUCTION 

In  all  of  our  previous  derivations  of  Snell’s  Law  of  Refraction,  our  functionals 
are  of  the  form: 

Ay]  =  [  +ya  dx . 

J  a 

This  represents  the  integral  of  a  function,  f(x,y)  with  respect  to  the  arc  length  s, 
where  ds  =  yj  1  +y/z  dx.  In  this  case,  the  Euler-Lagrange  equation  can  be  reduced  to 


the  following: 


A-//- 


where  f  =  f(x,y).  An  explanation  of  this  derivation  is  found  in  Gelfand  and  Fomin 
[Ref.  6:p.  19].  In  our  various  developments  of  Snell’s  Law,  the  function  f(x,y), 


represents  the  reciprocal  of  the  velocity  of  the  light  particles  in  question.  Recalling 
from  physics  that: 


sin4>, 

sin4>2 


-  =  nn  [Ref  1] 

V2 


we  now  introduce  the  term,  index  of  refraction:  nI2. 


Thus,  in  ail  of  our  earlier 


discussions  we  employ  a  functional  which  involves  the  index  of  refraction.  In  fact 
we  have  n(x,y)  =  f(x,y).  With  this  in  mind,  we  now  develop  certain  inverse 
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problems  to  illustrate  the  importance  of  the  index  of  refraction  to  various  acoustical 
and  optical  problems. 


B.  INVERSE  PROBLEMS 

We  pose  the  following  problem: 

Given  a  certain  host  medium,  find  the  index  of  refraction  which  will  allow  a  particle 
of  light  to  travel  from  a  point  A  at  (0,0)  to  a  point  B  at  (1,1)  or  to  a  point  C  at 

(y  - 1, 1)  along  the  following  curves:  (1)  a  parabola  (concave  down),  (2)  a  parabola 


(concave  up),  (3)  a  circle,  (4)  a  cycloid,  and  (5)  a  straight  line.  See  Figure  5.1. 


1 .  Parabola 

We  solve  these  different  problems  in  order,  starting  with  the  parabolas. 
Since  the  parabola  must  have  the  following  endpoints  y(0)  =  0,  and  y(l)  =  1  we 
necessarily  need  the  parabola:  y  =  x\  Applying  Euler's  equation  discussed  above 
we  arrive  at: 


or. 


-nxy  = 


,  _  ny' 


1  +y‘ 


~nx{  2x)  = 


2  n 

1  +4x2 


- n ,  _  1 

n  x(l  +4xz) 

which,  by  partial  fraction  decomposition  becomes: 
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Figure  5.1  Particle  Paths  in  an  Isotropic  Medium 
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This  integrates  to  yield: 


-*x  _  1  4x 
n  x  l  +4*2 


ln|n|  =  -In  |*|  +-^ln|l  +4*2|  +lnjC| . 


Therefore,  the  desired  index  of  refraction,  is: 


n(x)  = 


_  C/1T4I2 


The  other  parabola  (concave  up)  also  requires  the  same  endpoint  restrictions, 
y(0)  =  0  and  y(l)  =  1.  Therefore,  we  force  the  equation  of  our  parabola  to  be 

y  =  \[x.  Applying  Euler’s  equation,  we  are  left  with: 


-«/  =  -sl 

1  *ya 


X2fx~  l.± 


2  fx  yjx  (4*  -  1) 


from  which  it  follows  that: 


n  4*  + 1 


which,  after  integration  becomes: 


In  |n|  =  —In  |4*  + 1  j  +  ln  jC] . 

2 

So,  the  desired  index  of  refraction  for  this  parabola  is  the  following: 
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n(x)  =  Csj4x  +  1 . 


2.  Circle 

Now,  we  consider  the  circular  path.  Again,  we  desire  that  our  curve  pass 
through  the  points  A  and  B.  Therefore,  our  equation  for  the  required  circle  is: 

x2  +  (y-1)2  =  1 

which  yields: 


y  =  yj  1  -x2  +  l . 


From  our  earlier  discussion,  the  Euler  equation  is: 


from  which  we  get: 


1  «-y 


n\ 


I  -X2  1 


vr^"2 


1  +■ 


1-x2 


or. 


from  which  follows: 


x  -n(l  - x 2) 


or. 


-n 


\j  1  - x 1  \j\  - x 2 


Therefore,  after  integrating  we  arrive  at: 


ln|n|  =  -In  |*|  +ln  |C|  =  ln|— |  +ln  jCj . 

x 

Hence,  the  desired  index  of  refraction  for  this  medium,  in  order  to  travel  along  a 
circle  from  point  A  to  point  B  is: 

n(x)  =  . 

x 


3.  Cycloid 

Recalling  from  differential  equations,  the  parametric  equations  of  a 
cycloid  are  of  the  form: 


x  =  y(4>-sin<|>) 
y  =  y(l  -cos<t>)  [/&?/.  13]. 

Since  we  desire  that  our  cycloid  pass  through  the  points  A  and  C  we  necessarily 
reduce  the  above  equations  to: 

x  =  q>-sind> 
y  =  l-cos<J>. 


By  trial  and  error  we  determine  that  the  function  associated  with  a  cycloid  is  only 
a  function  of  y.  As  such,  the  corresponding  Euler  equation  is: 


ny" 

n’~° 


ny  1  2 y" 
n  2  [+/ 

The  chain  rule  and  integration  yield: 


(C!Uliy,(C^Ldy-.^^L&dx 

Tl  1  +yQ  1  +V72 


ln!«!  =  ln|  1  +  y^j  +  ln|Cj 


niy)  =  cjl+y* . 


Now  since 


*L  1  -cos<j> 


we  have: 


,  ,  1  +  sin24> 

niy)  -  C  | - ^ 

\  (1  -cos<|>)2 


niy)  =  C 


2(1  -cos<J>) 


\  (1  -cos<|))2 

Therefore,  the  desired  index  of  refraction  for  the  cvcloid  is: 


niy)  =  C .  - 

N  y 
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4.  Straight  Line 

Our  final  problem  involves  the  path  of  a  straight  line  from  point  A  to 
point  B.  That  is,  y  =  x.  Using  the  Euler  equation,  we  arrive  at: 


ny~n^'  ~~~H  =0 
1  +y 

/  ny 
-nv  =  — 

1  +ya 
n,  =  0 

which  integrates  to  yield  the  following: 

n(x)  =  C. 

Therefore,  the  index  of  refraction  in  the  case  of  a  straight  line  path  is: 

n(x)  =  C. 

Computing  the  various  indices  of  refraction  does  indeed  give  us  the 
desired  curve  from  a  point  A  to  a  point  B.  For  sake  of  brevity,  we  consider  the  first 
case,  that  of  the  parabola  (concave  down).  That  is,  we  show  that: 


iin./[>’]  =  f 1  *—X—  \)\  +ya dx 

Jo  x 


will  yield  the  curve  y(x)  =  x2. 

Applying  the  usual  Euler-Lagrange  equation,  that  is: 


Fy-J=-Fy,  =  0,  yields: 
y  dx  y 


d  Vi  +4x2  y' 

*.  *  {iTf 


58 


from  which  it  follow's  that: 


l/l  +4x2 

X 

y 12 

cV 

1  +ya 

1  +4x2 

„2„2 

v/2  -  _ 

C  X 

y - 

1  + 

(4-cV 

Setting  C  =  2  yields: 

v"  =  2x 

• t 

or, 

y( x )  =  x:  which  is  the  equation  for  the  parabola  (concave  down)  passing  from  point 
A  to  point  B  in  our  problem. 

C.  GENERAL  EXPRESSIONS  FOR  THE  INVERSE  PROBLEM 
1.  Function  of  X  Only 

It  is  possible  to  derive  a  single,  generic  expression  for  a  family  of  curves 
in  terms  of  any  function.  f(x,v),  where  the  various  functions  f(x,y)  determine  the  type 
of  curve  generated.  We  now'  show  the  development  of  such  an  expression.  We 
begin  with  our  standard  form  of  our  functional: 


J\y]  =  +ya  dx . 

We  first  restrict  our  attention  to  those  functions,  f  of  the  form  f  =  f(x),  that  is,  those 
functions  that  are  only  dependent  upon  x.  This  leaves  us  with: 

J[y\  =  fabf(x)JT^dx . 

Now,  since  the  integrand  does  not  depend  on  y,  the  corresponding  Euler  equation 
is: 

F,  =  C 
or, 

A*)y'  _c 

/TT? 

which  further  reduces  to: 


Mx)2-C2 


Now  integrating  the  expression  yields: 


-/ 


cdx 


\lfix)2-c2 


This,  then  is  the  simple  expression  we  have  been  looking  for  when  the  function  f(  x,y) 
is  only  dependent  upon  x.  With  this  expression,  we  can  now  input  any  function  and 
in  turn  generate  a  family  of  curves.  As  an  example,  we  make  use  of  an  earlier 
problem,  that  of  the  path  of  a  circle.  In  that  development,  we  saw  that  the 
corresponding  function  for  a  circular  path  was: 


fa)  =  - 

X 


For  simplicity,  we  will  let  C  =  1.  That  is,  our  problem  involves  flx)  =  —  and  we 

x 


wish  to  show  that  this  will  generate  a  family  of  circles. 


Starting  with  the  simple  expression  we  have: 

r  cdx 

y  =  — ■  — 
{jitf-c1 


which  becomes: 


Letting  u  =  1-C:x:  we  have: 


which  becomes: 


J  -  -c2 
N 


-cV 


>■  = — c  f— 

2  J  fu 


y~ci 


CJ 1  ~C2x 2 
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(y  -  C,)2  +x2  = 


_l_ 

c 2 


which  is  the  equation  of  a  circle,  centered  at  (0,C,)  with  radius 


This 


procedure  can  be  followed  for  any  function  and  will  generate  a  family  of  curves 
associated  with  that  particular  function. 


2.  Function  of  Y  Only 

We  now  turn  our  attention  to  those  functions  which  are  dependent  upon 
y  alone.  We  derive  an  expression  which  will  also  produce  a  family  of  curves.  As 
before,  we  start  with  the  now  familiar  functional: 


Jly]  -[tA^J>)\lUysdx. 

J  a 

We  consider  only  those  functions  dependent  on  y: 


J[y]  =  +ya  dx. 

J  a 

Since  the  integrand  has  no  explicit  dependence  on  x,  we  employ  the  following  Euler 
equation: 

F  -  y'Fy  =  C. 


So,  we  have: 


/tv)'/! 

or. 
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which  reduces  to: 


{l V 


or, 


y'  = 


N 


■Ay)2-c2 
c2 


After  integration  this  becomes: 


jt  + 


-/■ 


cdy 


v(/b02-' 


This  then  is  the  expression  which  produces  any  family  of  curves  dependent  only 
upon  the  function  f(y).  We  now  show  that  this  expression,  when  given  the  proper 
function  f(y),  yields  the  family  of  curves  for  a  cycloid.  Recalling  the  earlier 

rJy  1 

discussion,  the  corresponding  function  is  fly)  =  — — .  For  simplicity  we  let  C  -  —  . 


fy 


f2 


This  reduces  the  function  to  fly)  =  — .  Plugging  this  into  the  earlier  derived 


expression,  we  have: 


x  + 


««M 

J  rr 


cdy 


/r 


This  leaves  us  with  a  Volterra  integral  equation  of  the  first  kind: 
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x  + 


a,cf’-*L 

Jo  n — 

r,-' 


For  a  further  explanation  and  development  of  the  Volterra  integral 
Keener  [Ref.  25:p.  101].  Using  the  Fundamental  Theorem  of 
differentiating  both  sides  with  respect  to  y,  we  are  left  with: 


equation,  see 
Calculus  and 


or. 


dx 

dy 


This  reduces  to: 


or, 


dy  _  1 
dx  ~  ~C\ 


y-C2 


1 

C2)’ 


y(ya  + 


This  differential  equation  represents  the  equation  of  a  family  of  cycloids  [Ref.  27]. 


We  include  a  solution  to  this  differential  equation  for  completeness. 


Let  —  =  k2 
C 2 


and  let  y  =  k2sin:t.  Then,  we  are  left  with  the  following  expression: 

(1 +y'J)k2sin:t  =  k2 
or, 
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This  reduces  to: 


1  =  _L 


y'  =  cot(t).  We  consider  only  the  positive  square  root  in  this  case,  since  the  slope 
of  the  graph  is  positive  here.  This  then  yields: 

dx  =  —  =  tan(r)*fy . 
cot(r) 

However,  dy  =  k2(2sin(t)cos(t))dt.  Thus, 

dx  =  2k2sin2tdt. 

We  now  let  <l>  =  2t,  apply  the  initial  conditions,  and  make  use  of  the  double  angle 
formula,  sin2r  =  - — ^^,to  arrive  at: 


This  yields: 


dx  =  —  (1  -cos<|))  dfy 


dy  =  —  (sin<|>)  d< j> 


x  =  —  (<t>  - sin<J>)  +C, 


y  =  —  (l  -cos4>)  +  cz. 

Using  the  fact  that  at  =  0,  x  =  0  and  y  =  0  yields: 
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X  = 


*1 

2 


(4>  -sin<J>) 


and. 


y  =  y  (1  -cos4>) . 

These  are  the  equations  of  a  family  of  cycloids.  Again,  this  procedure  can  be 
followed  for  any  function  f(y)  and  will  yield  corresponding  families  of  curves. 
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VI.  CONCLUSIONS 


This  thesis  presents  several  derivations  of  Snell’s  Law  of  Refraction  and  one 
of  the  principle  of  reflection.  In  all  of  these,  we  place  different  constraints  upon  the 
problem  and  produce  the  familiar  laws  governing  reflection  and  refraction.  The 
main  result  of  this  research  is  the  expression  found  for  Snell’s  Law  of  Refraction  in 
which  the  velocity  of  the  light  particles  in  the  medium  depends  upon  direction.  We 
also  use  inverse  problems  to  find  the  index  of  refraction  needed  to  move  along  any 
desired  family  of  curves.  Finally,  we  derive  a  general  expression  that  generates  a 
family  of  curves  given  any  index  of  refraction. 

This  research  has  several  applications  in  the  area  of  acoustical  and  optical  wave 
propagation.  In  the  area  of  acoustics,  the  inverse  problems  enable  us  to  predict  the 
path  followed  by  sound  waves  for  a  given  medium.  This  would  allow  for  military 
use  in  the  area  of  prosecuting  submarine  contacts  and  in  ocean  bottom  surveys.  The 
inverse  problems  also  have  applications  to  the  bending  of  light  rays  used  in  laser 
technology. 

Although  this  paper  derives  several  generalizations  of  Snell's  Law  of  Refraction 
there  are  still  areas  that  require  further  research.  For  instance,  consider  the  inverse 
problem  in  space  where  the  index  of  refraction  is  a  function  of  three  variables  rather 
than  two.  Also,  in  all  of  the  derivations  in  this  paper  the  Hamilton-Jacobi  equations 
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were  not  used.  These  equations,  from  the  calculus  of  variations,  could  be  used  to 
derive  a  generalization  of  Snell’s  Law  of  Refraction. 
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